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Abstract 

We extend the polaron variational treatment previously developed for the propagator to the case 
where one nucleon and n external mesons are present. Using the particle representation of the 
scalar Wick-Cutkosky model this is done in lowest order of an expansion of the exact action 
around a retarded quadratic trial action. In particular, we evaluate the form factor for scattering 
of mesons from the scalar nucleon and determine the radius of the dressed particle. After analytic 
continuation to Minkowski space we study elastic meson-nucleon scattering both analytically and 
numerically near threshold and show that it is essential to incorporate the correct behaviour of the 
retardation function at large proper times. Only if this is done the optical theorem is approximately 
fulfilled over a range of energies and coupling constants. 



1 Introduction 



This is the third paper in a series which explores the use of variational principles in the particle 
representation of field theory. Following Feynman's treatment of the polaron [1] we have applied 
similar techniques [2, 3] to the simplest scalar field theory which involves a Yukawa interaction of 
heavy particles ('nucleons') with light mesons. This is the Wick-Cutkosky model [4, 5] described by 
the following Lagrangian in euclidean space time 

L = \ (d^f + \ M l& + \ (d^f + |mV " 9&<P W 

Here M is the bare mass of the heavy particle $, m is the mass of the light particle ip and g is 
the (dimensionfull) coupling constant of the Yukawa interaction between the two particles. In the 
present work we will concentrate on the Green function for one nucleon interacting with an arbitrary 
number of mesons. After integrating out the mesons (which is possible if the nucleons are quenched) 
we obtained in (I, II) the following generating functional for the corresponding (connected) Green 
functions 



Z[j,x] = const. / d(3 exp --M 2 / Vx{t) exp (-Sgff [x(t)] - S 2 [x(t), j] 

JO \ I J Jx(0) = 



/x(0) = 

Here the effective action for the nucleon is given by 



rP 1 n 2 rf j rP r rl A n 1 

and the meson source terms are contained in 

f (P ( d 4 a e iq '(.y- x ( T )) 

s^(r),j] = - g J^yml trj^^-r . (4) 

Eqs. (2 - 4) define the "particle representation" of the Wick-Cutkosky model in the 'proper-time gauge' 
[7] for the sector of the theory which we consider here: they are formulated in terms of trajectories 
x(t) of the nucleon which are parametrized by the proper time r and obey the boundary conditions 
x(0) = and x((i) = x. To obtain the Green functions for one nucleon (propagating from to x) 
and n external mesons one has to do the usual differentiation with respect to the meson sources j(y), 
perform the remaining path integral over all trajectories of the nucleon and finally integrate over 
(3 from zero to infinity with the weight exp(— (3Mq/2). It is, of course, impossible to perform this 
path integral exactly. In (I, II) we have therefore approximated it variationally on the pole of the 2- 
point function by a retarded quadratic two-time action. Various parameterizations for the retardation 
function which enters this trial action have been investigated and for the most general case we have 
solved the variational equations numerically. This was possible for values of the dimensionless coupling 
constant 

where M is the physical mass of the nucleon. For a greater than the critical value a c = 0.815 no 
real solutions could be obtained. The emergence of a critical coupling constant in our variational 
calculation points to the well-known instability present in the Wick-Cutkosky model [6]. 

It is the purpose of the present paper to extend the treatment given in (I, II) for the 2-point 
function to processes which involve n external mesons. We will do this in 'zeroth' order in the 
difference between the exact effective action (3) and the trial action 

AS = S e ff - S t , (6) 
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i.e. we will utilize the trial action S t which has been determined previously for the 2-point function 
to evaluate the (2 + ra)-point function. We will demonstrate that already this lowest order gives 
reasonable and non-trivial results. In particular, we will study in detail meson-nucleon scattering in 
this model which requires an analytic continuation of our results to Minkowski space. Remarkably 
this was already anticipated by K. Mano [8] who first pioneered the use of polaron variational methods 
in the Wick-Cutkosky model. 

This paper is organized as follows: in Section 2 we recall some elements of the polaron variational 
approach which we will need in the present paper. In Section 3 we consider the case where there are n 
external mesons and derive the lowest order approximation to the corresponding (2 + ra)-functions in 
an expansion of the exact effective action around the trial action. Section 4 discusses the special case 
of G 2) 2, the 'Compton' amplitude, the analytic continuation and the correct form of the trial action. 
In Section 5 we present our numerical results for differential and total cross sections near threshold. 
The main results of this work are summarized in the last Section, whereas some technical details are 
relegated to the Appendix. 



2 Polaron Variational Approach 

Following Feynman's treatment of the polaron problem we have performed in (I) a variational calcu- 
lation of the 2-point function with the quadratic trial action 

rP 1 f 13 f Tl 2 

S t [x] = / dr -x 2 + dr x \ dT 2 f{T 1 -T 2 ) |>(ri) - x(t 2 )] . (7) 
Jo Z Jo Jo 

Here /(ti — r 2 ) is an undetermined 'retardation function' which takes into account the (proper) time 
lapse occurring when mesons are emitted and absorbed on the nucleon. In actual calculations it is 
more convenient to use the Fourier space form 

oo 

S t [b]=J2A k b 2 k , (8) 

k = 

where the b k are the Fourier components of the path x(t) and the Fourier coefficients A k are either 
parametrized or considered as variational parameters. Eq. (8) contains as a special case the free 
action for which all A k equal unity. Near the nucleon pole the 2-point function should behave like 1 

where < Z < 1 is the residue. As was shown in (I) this requires the proper time (3 to tend to 
infinity. For (3 — > oo all discrete sums over Fourier modes A k turn into integrals over the 'profile 
function' A(E = kir/(3). Mathematically this is achieved by using Poisson's summation formula [9] 
which, for an even function F(kir/(3), reads 

00 ( kir\ Pi r°° 1 9/? 00 r°° 

£ F hd=-/ dEF(E)--F(0) + ^ £ / dEF(E) cos(2n/3£) . (10) 
t^i V P J k Jo 2 ^ Jo 

For (3 — > oo the terms in the last sum are exponentially suppressed. 

The profile function A(E) is linked to the retardation function f(cr) through 

8 f°° Err 
A(E)= l + ^l daf(a) sin 2 — . (11) 



'In this work all 4-momenta are taken to be euclidean. 
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In Refs. (I, II) we have studied the 'Feynman' parametrization given by 



w 



w 



(12) 



which leads to 



and an 'improved' parametrization 



A F (E) 



v 2 + E 2 
w 2 + E 2 



w 2 1 



which gives 



AjiE) =1 + 2 



w 



wE 



2w 



E 

arctan — 
w 



w 
2E 



E 2 

1 + -7 



(13) 



(14) 



(15) 



In both cases v, w are parameters whose values have to be determined by minimizing the variational 
functional given in Eq. (27) of (II). As the bare mass M requires renormalization and the physical 
mass M = 939 MeV is fixed by the pole position of the propagator, the minimum does not have a 
physical meaning but only tells us how good the variational ansatz is. The numerical calculations in 
(II) have shown that the 'improved' parametrization is significantly better than Feynman's because it 
includes the correct short-time behaviour. Note also that the singularity structure of the corresponding 
profile functions in the complex E'-plane is quite different: whereas A F (E) exhibits poles at E = ±iw, 
the profile function of the 'improved' parametrization has branch points at E = ±iw. 

Within the gaussian ansatz the best approximation is obtained by not imposing any specific form 
for the retardation function but determining it from varying the variational functional with respect 
to the profile function A(E). As shown in (I) this gives the following expression for the 'variational' 
retardation function 



/var(^) 



1 



32tt 2 fi 4 (a) Jo 



du 



m 2 „, ,1-u X 2 M 2 a 2 
1 H u (<r) r-— - 



( . . XMa \ . . 

e [mfi{a),-^—^,u j . (16) 



Here < A < 1 is a variational parameter determined by solving the coupled system of variational 
equations. In addition, we have used the abbreviations 



e(a, b, u) 



a 2 1 



exp ( - - 



u b 2 
— — — u 
2 



(17) 



and 



, 4 f°° , 1 sin 2 (£W2) 



E 2 



:is) 



We call j-i 2 {cf) the 'pseudotime' because 



(19) 



and 



T 



(20) 



where A = A(0) is the value of the profile function at E = 0. It should be noted that this 
variational retardation function has the same l/<7 2 -behaviour for small relative times as the 'improved' 
parametrization (14). 
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3 Green Functions with n External Mesons 



In the previous work we have been considering the propagator of the nucleon dressed by (virtual) 
mesons and have used this Green function to fix the variational parameters appearing in the trial 
action. It is of course possible to apply the same procedure for any other Green function as well. 
Alternatively, one can estimate other Green functions by evaluating them with the variational pa- 
rameters determined from the nucleon propagator. We shall follow this procedure for the general 
Green function with two external nucleon and n meson legs which we shall call the (2 + ra)-point 
function. With the help of the reduction formulas these Green functions describe the matrix elements 
for processes like meson absorption on a nucleon, meson-nucleon scattering or nucleon-antinucleon 
annihilation in our simplified scalar model. Nonperturbative production of multiboson states has 
attracted a lot of theoretical interest recently (see e.g. [10]), in a different context. 

For simplicity we will work to zeroth order in (AS) and within the framework of the 'coordinate 
averaging' described in (I). The general expression for the (2 + ra)-point function in coordinate space 
may be obtained from the generating functional for connected Green functions (Eq. (2)) by repeated 
differentiation with respect to the source j(y). After Fourier transformation we may write 



G 2 , n (p,p'; {q}) 



const. / dftexp 



2 



x(fi)=x 



Vx(t 



x(0) = 



n 



d 4 xe %p '* x 



g / d Tl e^< T ' 



s eff Kr)] 



(21) 



where the external meson propagators as well as an overall momentum conserving delta function have 
been removed, while the truncation of the external nucleon propagators still remains to be performed. 
Our convention is that nucleon lines are ingoing with momenta p and p' whereas mesons are outgoing 
with momenta g, . The effective action S e ff for the nucleon is given in Eq. (3). We shall determine 
the normalization of Eq. (21) later. The proper time r, appearing in Eq. (21) may be interpreted 
as the time at which the meson with momentum g, couples to the nucleon, the latter being at the 
space-time point a;(r,). 

To zeroth order in (AS), G^^iPiP 1 ] {q}) is given by Eq. (21), with S e ff replaced by S t . Using the 
definition for S t in Eq. (8) one may do the integral over the endpoint 



d 4 x exp 



-ip' ■ x - A ^ + 



as well as the path integral 

r x(p)=x 



Vx(t) exp 



a;(0) = 



where 



CO 

(iK ■ h ~ A k b 

2^ 





'27T/3' 


2 


1 




[ A \ 


exp 


2A (3 



\qiTi -p'(3 



k = l 



const. 

(27T/3) 2 



exp 



oo \f2 



E 

k = l 



4A k 



k = l A k 



At 



irk 



n t 
^-^ KTTTi 

2^ sin — — qi 



(22) 



(23) 



(24) 



One may perform the sum appearing on the right hand side of Eq. (23) using the summation formula 
(10). One obtains, up to exponentially small terms in (3, 

f n \ 2 n 

E Tiqi + VJ7^ X) qi ' q i mil 



oo \ /2 
A h 



E 

k = l 



1 



2/3A(0) 



dE ( 1 



1 f 00 dE ( 1 1 \ 

+" E * • I {a{E) ~ 1(0) J Sin( ^ } Sin( ^ } • 



(25) 
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Figure 1: The definition of the relative times and momenta. Note that r = 0, r n+1 = (3, pi = — p and 

Pn + l = P'- 



Note that the first term, which is quadratic in the proper times r and therefore potentially troublesome, 
will cancel against the identical term appearing in Eq. (22). 

The expression (21) for the ra-point function involves an integral over all possible orderings of the 
times {t;}. Equivalently, we may relabel the integration variables such that they correspond to the 
particular ordering < t x < r 2 . . . < r n < (3 and sum over all possible permutations of the momenta 
{qi}. Furthermore, it is more sensible to use relative rather than absolute times, so we shall define 
(see Fig. 1) 



a 



Q!i 
n + l 



T 2 - Tl 
7"n 7"n — 1 



(26) 



so that 



n + l 



T~i = ^ Oi k and fi = ^ a k . 

k=l k=l 

The Jacobian of this transformation is one, and the product of integrals simplifies to 

I'fl i>T n i'T-2 n + l ^oo 

d/3 / dr n / dr n _i . . . dri = TT / da t . 

Jo Jo Jo i=1 Jo 

With this transformation the terms linear in the times appearing in the exponential become 



1 



2A 



[i p' 2 - 2p' ■ liTi + * ' 1i mm ( r i' T i) 



1 



n + l 



24 T, a >P* > 

ZA ° i=i 



(27) 



(28) 



(29) 



where the momenta p, are defined in Fig. 1 and we have used the fact that the overall momentum is 
conserved. We are thus led to the following expression for the (2 + ra)-point function: 

n + l 



G 2 ,n{p,v'; {q}) 



^ n-t-± „oo 

T~n ^ II \9 / da i eX P 



'2 An 



{p 2 + A M ' 



k = 



1 



6 



1 " r°° dE \ 1 1 



(30) 



• sin 



^ X] a < sin \ E X! a ' 



(=1 



(=1 



This expression is valid for all non-negative integers n and, as we shall see, is properly normalized. 
Before treating the case where n / 0, let us first consider what the result is for the nucleon propagator 
at this order. 

3.1 n = 

Applying the Poisson summation formula (10) to the infinite product of Fourier coefficients gives 



The remainder 



II T~2 = IT exp — / dE lo % A ( E ) + Ex (/ 3 ) 

fc=o A k A I k Jo J 



ME) 



sm{2n(3E) 



(31) 



(32) 



is exponentially suppressed for large (3. Hence close to the nucleon pole, i.e. for (3 — > oo, we obtain 



p 2 + A M 2 + 



4A n 



7T Jo 



G 2 fi{p,p) -- 

So at zeroth order the residue of the nucleon propagator is one 

ZM = 1 

and the physical mass is given by 



dE log A(E) 



M 2 



A 



Mo +- / dE \og A{E) 
vr Jo 



(33) 



(34) 



(35) 



Note that the singular behaviour of the 'improved' and the 'variational' retardation functions for small 
a leads to an 1/E fall-off of the corresponding profile function for large E. This in turn requires (even 
in zeroth order) an infinite renormalization of the bare mass M since the ^-integral in Eq. (35) does 
not converge. 

As we have already noted, Eq. (33) is only valid near the pole. It is usually not possible to give a 
closed expression for the infinite product in Eq. (31) or the infinite sum in Eq. (32) which is also valid 
for subasymptotic /3's. One case, however, where this may be carried out exactly is for the Feynman 
parameterization given in Eq. (13), in which case we have 



1 



k = 



w sinh wfi 
v sinh vfi 



(36) 



so that the propagator becomes 



£2,0 (P,p) 



w 



d(3 



sinh wfi 



sinh vf3 _ 



exp 



JL 

2A 



(p 2 + A Mo 2 



(37) 
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Figure 2: Singularity structure and integration contours for the 'improved' profile function in the 
complex E'-plane. 



One may even carry out this last integral over (3 to obtain 

,2 



G 2 ,o(p,p) 



w 
4 v 3 



2w 
a -\ 

v 



1 * 



2w 
a -\ 

v 



+ (a - 1) * (a) 



a H 1 

v 



w 

a -\ 



(38) 



with a = w 2 (p 2 + M 2 )/(4u 3 ) . Here is the digamma function [11], which has simple poles at 

a; = 0,-l,-2,.... 

It is clear that this Green function has the unfortunate feature of having not only the pole situated 
at the nucleon mass, but also an infinite sequence of 'ghost poles' (unless v = w, which corresponds 
to g = 0, in which case all ghost poles cancel). Note that the residues at these ghost poles cancel each 
other, which is consistent with = 1 noted above. The closest pole is situated at p 2 -\-M 2 = — Av 2 /w, 
which numerically turns out to be more than about (2M) 2 away from the nucleon pole. Hence, for 
the purpose of describing low-energy processes, these ghosts seem not to be relevant physically. 

What happens to these ghost poles for a general profile function A(E) ? Without a specific 
parametrization it is possible to investigate the nearest singularities of the zeroth order propagator 
by keeping only the n = 1 term in the sum for the remainder (32) and expanding the exponential. 
This gives 



£2,0 (P,p) 



1 



p 2 + M 2 



+ 



8A 



dEE 



A'(E) 



1 



(39) 



A(E) (p 2 + M 2 ) 2 + (4A E) 2 ' 

To make further progress we need some knowledge of the analytic behaviour of A(E) in the complex 
E'-plane. It has already been observed that the 'improved' parametrization of the profile function has 
branch cuts on the imaginary axis starting at ±iw. Assuming this structure to hold in general (with 
the cuts starting at some value HE ) we can, after making use of the fact that the profile function 
is even in E, deform the integration contour as shown in Fig. 2 and obtain 

r(E) 



G 2 ,o(p,p) 



1 

p 2 + M 2 



+ 



dE 



p 2 + M 2 + 4A E 



(40) 
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with 

V ' 7T A{iE+€) ^ ' 

Eq. (40) has the form of a Kallen-Lehmann representation [12] for the euclidean propagator and shows 
that the ghost poles in Feynman parametrization turn into branch points of the off-shell propagator 
- at least for the 'improved' parametrization. It is tempting to identify these branch points with 
the beginning of the meson production cuts which should be present in the exact 2-point function. 
However, the threshold at p 2 + M 2 = — 4A E turns out to be too high if the numerical values of 
A and E = w tabulated in (II) are inserted, and, more seriously, the weight function r(E) is not 
positive definite. The latter deficiency is also obvious from the result 

[CO 1 [CO J 

/ dEr(E) = — / dE— [In A(iE + e) - In A(iE- e) ] = 0. (42) 
Je Je dE 

Although this is fully consistent with the general sum rule for the weight function in the Kallen- 
Lehmann representation [12] and = 1 , it means that, in general, we have 'ghost branch points' 
instead of ghost poles. These unwanted properties are not as disastreous as it seems at first sight: 
it should be remembered that we have extrapolated away from the nucleon pole with the variational 
parameters fixed to their on-shell value. As shown in (II) the off-shell variational equations necessarily 
lead to a dependence of the variational parameters on the virtuality p 2 + M 2 which will cure these 
deficiencies. We will not pursue this in the present work but will keep the variational parameters as 
determined on the pole of the 2-point function. 



3.2 n ^ 

Having determined the residue of the nucleon propagator at zeroth order, we may now check that the 
Green function (30) is properly normalized by considering the small coupling limit. In order to obtain 
the lowest order contribution in the coupling g, we just need to set the profile function A(E) equal to 
one and M to M. One can now do the integrals over the times a, to obtain 



c£ e< W;{?}) 



i 

~9 



E 



n + 1 

n 



2g 



V({q}) i = l 



P % 2 + M 2 



(43) 



as required. 

Let us now specialize to the case where the nucleon is on-shell. As discussed in (I) this corre- 
sponds to considering the integration region where the corresponding time differences oix and a n+ i 
go to infinity. Hence, using again the Poisson summation formula, the terms quadratic in the meson 
momenta may be simplified by writing the product of the two sine functions as a sum over two cosines, 
one of which only gives rise to exponentially small (in the variable o^) terms. One obtains 



1 



where 



dE 



1 



1 



7T Jo E 2 IA(E) A(0) 



sin 



e ^2 a i sm \ E ^2 a ' 



i=i 



i=i 



l r°° dE 
2^ Jo ~EP 



(max(ij) 
;=min(ij') + l . 



cos(£'r) 



(44) 



(45) 



-ME) A(0). 

Note that the argument of £ in Eq. (44) is zero if i = j. Furthermore, r does not depend on oix or 
cn n+ i for any values of i and j and so we may now carry out the integrations over these two times, 
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(a) 




(b) 



(c) 



Figure 3: Some diagrams included in the (2 + n)-point function. Each solid line stands for a fully 
dressed nucleon propagator. 



yielding the propagators for the external nucleon legs. The (2 + ra)-point Green function for on-shell 
nucleons may therefore be written as 



where 



G 2 , n (p,p'; {q}) 



p 2 + M 2 p,2 + M 



G? n (p,p';{q}) , 



(46) 



r poo 

G 2 tT n(PiP''Aq}) = 2 ^o ^ n(«| dai exp 



V({q})i = 2 



~2A, 



{n I max(ij') 

i,i = l \;=min(ij') + l 



(47) 



and we have used the zeroth order expression Eq. (35) for the physical mass. 

Note that all momenta in Eq. (47) appear quadratically in the exponent. This is due to the fact 
that it corresponds to the zeroth order term in a variational calculation with a quadratic trial action. 
Nevertheless, already at this order Eq. (47) contains a remarkable amount of information. To see this, 
it is instructive to examine its behaviour in perturbation theory. It suffices to expand the exponential 
in the g,'s. The leading term corresponds precisely to the tree diagram already encountered before 
(Eq. (43)), with fully dressed nucleon propagators for the internal lines and bare vertices for the 
external pions. An example of this type of term is shown in Fig. 3 a. Next, there is a class of 
diagrams containing only diagonal terms (i = j). These correspond to similar diagrams, however this 
time with one or more dressed vertices (e.g. Fig. 3 b). This class of diagrams is reducible on each side 
of each dressed vertex. This is due to the fact that £(0) is independent of the a's and so the integrals 
over the times may be carried out, yielding a product of propagators and form factors. Furthermore, 
there is a class of diagrams where off-diagonal terms (i / j) are also present. Because for these terms 
£ contains at least one a, these correspond to diagrams which are not reducible on each internal leg 
(e.g. Fig. 3 c). For example the term in the expansion with i = 1 and J = 2 gives qi ■ q 2 £(0^2) an d 
thus corresponds to a diagram with meson exchange across the first and second vertex. Lastly, note 
that the exponential gives rise to any combination of the above, to all orders of the coupling. 

The function £(r) which we have introduced in Eq. (45) is closely related to the pseudotime /U 2 (t). 
Indeed using Eq. (18) we obtain 



= ?(0) + 



4A(0) 



(48) 
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(a) (b) 

Figure 4: Third-order graphs for Gr 2 ,i : (a) self-energies, (b) vertex correction. 



Note, in particular, that 



A, 



1 + 4£'(0) 



(49) 



so that all on-shell Green functions are completely determined by the function £(r). For the Feynman 
parametrization (13) we have the analytical result 



Mr) 



v - w 



4v 3 



(50) 



which illustrates the fact that £(r) is an exponentially decaying function of r (see Eq. (A. 27)). 

Finally, for convenience we explicitly write down the lowest order Green functions, where all of 
the above points may be illustrated : 



G , 2 t 2(p,p';gi,g 2 ) 



2gA e-« 2 «°) 
2g 2 A e-M + &W 
a 2 



da 2 



exp 



^-(p 2 2 + M 2 ) - 2q x ■ q 2 £(a 2 )] + (ft f+ q^j 



(51) 



(52) 



G , 2 t 3(p,p';gi,g2,g 3 ) 



2g 3 A e-M +qa > + &W V / 



dot-) da-* 



?{{<!})■ 



" 2 -(P2 2 + M 2 ) -^(p 3 2 + M 2 ) 



6XP L 2A ^ ' ' 2A 
• exp [-2gi • q 2 £(a 2 ) - 2q 2 ■ q 3 £(a 3 ) - 2q x ■ q 3 £(a 2 + a 3 )] 

The expression for the (2 + l)-point function in Eq. (52) may be written as 

G^p'-q) = 2g e $F(q 2 ) , 



(53) 



(54) 



showing that the effective coupling constant for the meson-nucleon vertex is enhanced due to vertex 
corrections : 

#eff = 9 A > g . (55) 

This is qualitatively (but not quantitatively) similar to what one obtains from expanding the one-loop 
effective potential V^^, <p) (see II) up to the power of Q 2 (p . 
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Not surprisingly the elastic meson-nucleon form factor 

F(q 2 ) = e~^ 2 (56) 
is gaussian and defines the mean square radius of the dressed particle to be 

^> = 6 ^ = U? E i>(m-m)- < 571 

For the Feynman parametrization we obtain from Eq. (50) 

3 v 2 — w 2 



„2\ 



r 



F 



2 v 3 



(58) 



which is a well-known result in the polaron literature [13]. Table 1 gives the numerical values obtained 
with the different parametrizations for the profile function. We also have included the perturbative 
result from the vertex correction in Fig. 4 b which reads 

<r 2 >pert = {~ 2 I du [{1 _ u)m2 + um2 f ( 59 ) 



Table 1: Root-mean-square radius (in fm) of the dressed particle from Eq. (57). 'Feynman' signifies 
the result in the Feynman parametrization whereas 'improved' refers to the improved parametrization 
from Eq. (15). The radius calculated with the solution of the variational equations is denoted by 
'variational'. For comparison the perturbative result is also shown. 



a 


'Feynman' 


'improved' 


'variational' 


perturbative 


0.1 


0.014 


0.026 


0.018 


0.047 


0.2 


0.021 


0.038 


0.027 


0.066 


0.3 


0.028 


0.050 


0.034 


0.081 


0.4 


0.034 


0.062 


0.042 


0.094 


0.5 


0.041 


0.073 


0.050 


0.105 


0.6 


0.049 


0.088 


0.060 


0.115 


0.7 


0.060 


0.108 


0.073 


0.124 


0.8 


0.081 


0.149 


0.098 


0.133 



The values of the root-mean-square radius (Table 1) for the variational calculations differ much 
more from each other, as they do from the perturbative results, than the corresponding results for 
the residue in (II), even for small coupling. The reason for this difference at small coupling is that 
we have only calculated the (2 + 1)- (and higher) point function at zeroth order in (AS). So, by 
construction, they are only 'forced' to agree at the tree level in perturbation theory (which they do). 
If one calculates them to first order in (AS) the agreement for small couplings will be similar to the 
agreement for the residue of the two-point function determined in (II). Still, even at zeroth order the 
root-mean-square radii are quite consistent with each other. It should be noted that the radius is 
extremely sensitive to the behaviour of the profile function at small E, so the differences between the 
radii is a reflection of the differences between the profile functions seen in this region. 
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The radius (57) is due to the recoil of the nucleon when it absorbs the external meson while other 
mesons have been emitted into the meson 'cloud' surrounding the dressed particle. Since there is 
no direct meson-meson interaction the natural length scale for the radius therefore is the Compton 
wavelength 1/M = 0.21 fm of the nucleon. The actual values are, of course, dependent on the strength 
of the coupling and, as seen in Table 1 the root-mean-square radius increases with a. In principle for 
some large enough coupling it should become comparable to the empirical charge radius of the proton 
( ~ 0.8 fm ), were it not for the instability which forces the coupling constant to be less than a c . This 
is probably an unreasonable comparison in any case - we are treating the (bare) nucleon as a point 
particle here, so the root-mean-square radius is only due to effects from the meson cloud. 



4 Meson-Nucleon Elastic Scattering 

We now turn to a more detailed discussion of the 'Compton' amplitude given by Eq. (53). It consists 
of two parts, the direct diagram and the crossed one shown in Figs. 5 a and 5 b where we have 
set q = —qi, q' = q 2 for the meson and p' — > —p' for the outgoing nucleon. In terms of the usual 
Minkowski space Mandelstam variables 



-(q + pf 



t=-{p-p'f 



(60) 



we have for the crossed diagram 



G 



crossed 



2,2 



(s,t) = 2g 2 A e 2 ™ 2 «°) / da 2 exp 



and for the direct diagram 



«2 

'2A 



{s + t-M 2 - 2m 2 ) - (2m 2 - t) f (a 2 ) 



(61) 



G 2 f ect (s,t) 



2g 2 A 



2m 2 {(0) 



da 



2 exp 



«2 

"2Ao 



{-s + M 2 ) - (2m 2 -t)^(a 2 ) 



(62) 



It is worthwhile to point out some similarities and differences of these formulae with the usual skeleton 
expansion of higher Green functions: after the substitution r = a 2 /(2A ) in the proper-time integrals 
we indeed can identify the enhancement factor in front of the integrals as G\ 1 (q 2 = —m 2 ), i.e. the 
square of the meson-nucleon vertex function (52). However, the zeroth order propagator G 2<0 studied 
previously does not appear in Eqs. (61, 62) and consequently, none of its off-shell deficiencies will 
show up. This is because we only consider truncated Green functions which, by definition, are on the 
mass shell with respect to the external particles. The variational principle then yields well-behaved 
expressions even though the internal particles may be off-shell in a diagrammatic expansion. 

When applying Eqs. (61, 62) to meson-nucleon scattering the only question which remains is the 
convergence of the integral representations for the 'Compton' amplitude. Whereas Eq. (61) for the 
crossed diagram converges for s + t > M 2 + 2m 2 , which includes the physical region, this is not the 
case for the direct diagram. Indeed, Eq. (62) is only defined for s < M 2 and must be analytically 
continued for the physical values s > (M + to) 2 . 



4.1 Analytic Continuation 

To perform this analytic continuation we must investigate the analytic properties of the function £,(a 2 ) 
which has been defined in Eq. (45) as the Fourier cosine transform of the inverse profile function. 
Using the fact that the profile function is even in E and assuming that it has the analytic structure 
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(b) 



Figure 5: 'Compton' amplitude for meson-nucleon scattering in the zeroth order variational calcula- 
tion: (a) direct diagram, (b) crossed diagram. 



shown in Fig. 2 we may again deform the integration contour and express ^(a 2 ) by an integral running 
along both sides of the cut 

f oo 

dEp(E)e- Ea2 (63) 



where 



e(« 2 ) 



P{E) 



1 



2ttE 2 



Im — — . (64) 

The explicit form of the weight function p(E) for the various parametrizations can be found in the 
Appendix. In deriving Eq. (63) we have tacitly assumed that the profile function A(E) has no zeroes 
in the upper plane which can be verified for the 'improved' parametrization. Of course, if this is not 
the case one has to include the contribution from the corresponding poles in the integrand for ^(o^) • 
The integral representation (63) can now be used to evaluate the a 2 -integral for the direct diagram. 
After expanding the exponential and integrating term by term we obtain 



G 



direct 



(s,t) = 2/A e 2m2 « 0) 



2A n 



M 2 



+ 



E 



(t - 2m 2 )' 



dE\ ... dE n 



E 



piE,) ... p{E n ) - 



-L- {M 2 -s) + E 1 + ...E n 



.(65) 



Although Eq. (65) has been derived for s < M 2 , it now can be used to define the direct amplitude 
for other values as well, in particular on the upper side of the s-cut ( s — > s + it ) which corresponds 
to the physical region. The infinite series (65) may then be resummed by means of the representation 



a — it 



i I (It e~ 

>o 



ir(a-if) 



This gives 



G 



direct 



2,2 



(s,t) = 2g 2 A e 



2m 2 {(0) 



(It exp 



(66) 



(67) 
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Of course, Eq. (67) could have been derived more easily from Eq. (62) by the simple replacement 

a 2 = it , (68) 

which amounts to the usual Wick rotation back to Minkowski proper time. However, the derivation 
makes clear the crucial role played by the analytic properties of the profile function A(E) which 
determines the function £(r). In addition, the expanded version (65) of the direct amplitude is 
very useful to study thresholds and the threshold behaviour of the imaginary part of the scattering 
amplitude. This will be done in the next subsection. 



4.2 Threshold Behaviour and Optical Theorem 

As the zeroth order residue is one, the truncated Green function equals the scattering amplitude 

T(s,t) = G$(s,t). (69) 

This may be verified by calculating the Born terms in perturbation theory. Using, for example, the 
conventions for the Feynman rules advocated by Muta [15] 2 one obtains 



T Born = {2gf 



+ 



M 2 + {p+q) 2 M 2 +{p-q l ) 2 _ 



(70) 



which is identical with Eq. (43). 

The differential cross section in the center-of-mass system can be expressed as (Ref. [14], p. 245) 

% = *k™\' (71) 

and the total cross section by means of the optical theorem as 

CTtot(s) = 2]^ ImT(M = 0) - (72) 

Here p is the center-of-mass momentum. The Born amplitude (70) is, of course, purely real thereby 
leading to a vanishing total cross section and a violation of the optical theorem. The direct amplitude 
in our zeroth order variational calculation, however, develops an imaginary part after the analytic 
continuation has been made. Indeed, using the expanded version (65) of the direct amplitude with 
s — > s + it we have (in the physical region) 



00 ( — O m 2 ) n r°° 

Im T(s, 0) = 2nc/A e 2 ™ 2 «°> £ ^ ^ / dE x ... dE n p^) ... p{E n 

n = l n - Je " 



' 6 {A' {M2 ~ S)+El+ - En ) ■ (73) 

The 5-function may be used to perform, say, the ^-integration which gives a nonzero result only 
if the singularity lies within the integration interval. Since all integrations over the E^s start at E 
where the branch point of the profile function is located, we find that the n - th term in Eq. (73) only 
makes a contribution to the imaginary part of the scattering amplitude if 

s > thresh = M ' 2 + 2ra ^o£o n = 1,2... (74) 



2 These conventions produce the minimal number of i's, do not change the sign of the propagator when a Wick 
rotation is made and directly yield the transition matrix element. Our Lagrangian (1) gives a factor 2g at each vertex. 
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In particular, between the elastic and the first inelastic threshold for meson-nucleon scattering only 
one term in the sum (73) contributes and we have 

"» n>, 0) = "4*9 V.4. «» *> ( ^ ) ■ ^sh < • < ^Lsh ' < 75 > 

We see that the cut structure of the profile function A(E) not only determines the threshold position 
(74) for meson-nucleon scattering but also directly the total cross section. 

Let us now discuss qualitatively the effects produced by the different parametrizations which we 
have used in the variational calculation on the pole of the 2-point function. We immediately recognize 
that the rather small deviations which these parametrizations have given for the self-energy lead to 
very large ones for the scattering cross section. For example, the Feynman parametrization has no cuts 
at all but only poles and zeroes in the profile function. Consequently there is only a very pathological 
(and unphysical) scattering cross section when this parametrization is continued to Minkowski space: 
from Eq. (A. 2) in the Appendix we see that the total cross section would be a sum of 5-functions ! 

The 'improved' parametrization fares better at first sight: it has a cut in the upper E'-plane 
starting at E = iE = iw . Therefore 

^thresh " M = + 2AoEo ~ M ( 76 ) 

which should be equal to the meson mass m. Using Table 2 in (II) the right-hand side of Eq. (76) 
is found to vary between 408 MeV at a = 0.1 and 275 MeV at a = 0.8 . Although this is a far cry 
from the correct value 140 MeV it could be easily cured by performing the minimization with the 
constraint that the correct threshold is reproduced. A more serious problem, however, is that near 
threshold the weight function p(E) behaves like (see Eqs. (A. 5) and (A. 10) ) 

p(E) Cj (E - w) (77) 

with a positive constant Cj. In this kinematical region the center-of-mass energy is 

. = ( ^WT^ + V^TV-)' p - 4L sh + «• ^S 1 ■ < 78 » 

Using the threshold condition (76) we thus find that the weight function in Eq. (75) behaves like 

/ s - M 2 \ p^o p 2 (M + to) 2 
9 V. 2A ) ~^ 1 2A Mm 

This has two disastrous consequences: firstly, the imaginary part of the amplitude vanishes like p 2 
(i.e. not like |p|) leading to a zero total cross section (72) at threshold and, secondly, away from 
threshold ct^qI becomes negative when calculated by means of the optical theorem ! 

We will see in the next subsection that both these failures are not due to the zeroth order variational 
approximation but are a consequence of an inadequate form of the retardation function. 

4.3 An Extended Parametrization of the Retardation Function 

It is not surprising that many of the parametrizations for the retardation function investigated in 
(I, II) turn out to be inadequate for describing meson nucleon scattering. By means of a simple 
approximation we have found in (I) that the variational principle for the self-energy is sensitive 
mostly to small values of the proper time difference a. On the other hand, a scattering amplitude 
near threshold probes large values of this variable as can be inferred, for example, from Eq. (67). 
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(79) 




Figure 6: Retardation functions f(cr) versus proper time a for various parametrizations: 'Fey' de- 
notes Feynman's parametrization, 'imp' the 'improved' parametrization and 'var' the solution of the 
variational equations. Note that for large a the variational retardation function becomes negative. 



Fortunately the variational principle tells us what the best form for the retardation function is for 
all values of a, including a — > oo. Indeed, performing an integration by parts in Eq. (16) one obtains 



/var(o") 



9 



1 



32tt 2 n A {a) 



exp 



A 2 MV 2 
2 M 2 (a) 



m 2/ N 



1 / XMa 
au e muicr), — — — , u 
V (J.{(t) 



(80) 



Replacing /2 2 (a)/a by some average value one sees that the 'improved' parametrization just corre- 
sponds to the first term in Eq. (80). However, the second term dominates at large a and makes the 
retardation function even negative. This is shown in Fig. 6 where the retardation functions investi- 
gated in (II) are plotted as function of a . It can also be demonstrated analytically in the following 
way: recalling the asymptotic limit of the pseudotime given in Eq. (20) we evaluate the M-integral by 
Laplace's method 



du e 



-E(u)a 



27T 



<tE"(u ) 



-E(u )e 



(81) 



Here u is the value which minimizes the function E{u) in the interval between and 1. For the case 
at hand (see Eq. (17) ) 



E(u) 



m 



2A 



1 \ A X 2 M 2 

1 + u 

u J 2 



so that 



E u 



m 



AnXM ) 



XMm 



m 



E"(u ) 



A 2 n X 3 M 3 



> 0. 



m 



(82) 



(83) 



iQ/\iK/ / 2Aq 

Since E(u ) is smaller than the value A 2 M 2 A /2 which governs the exponential fall-off of the first 
term, we see that the variational retardation function asymptotically behaves like 



/var(^) 




-.E(llo) cr 



(84) 
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Evidence for a sign change of the variational retardation function for large a can already be found in 
the 'dip' of Ayar(-E') near E = which has been observed in the plots given in (II). As a side remark we 
note that the threshold for meson-nucleon scattering determined from Eqs. (84) and (83) is now much 
closer to the physical value than with the 'improved' parametrization: with E = E(u ) and Table 3 
of Ref. (II) the right-hand side of Eq. (76) now varies between 120 MeV at a = 0.1 and 99 MeV at 
a = 0.8 . Far more important than the remaining discrepancies (which may be easily removed by a 
constraint calculation) is the sign of the asymptotic retardation function and the fractional power of 
a which accompanies the exponential decay. This is because a retardation function of the asymptotic 
form 



c_ 



-E a 



leads to the following behaviour of the weight function near the branch point E 

p(E) oc (E-Eoy- 1 , 



(85) 



(86) 



the proportionality constant having the same sign as the constant C in Eq. (85). This is demonstrated 
in the Appendix. Thus the asymptotic behaviour (84) "kills two birds with one stone" : it produces a 
positive total cross section and the correct threshold behaviour of the imaginary part of the scattering 
amplitude to make the total cross section finite at p = 0. In addition, it leads to a lower minimum 
of the variational functional thereby showing that such a form is a consequence of the variational 
principle applied to the present model and not an arbitrary parametrization. 

We incorporate the behaviour of the variational retardation function at small and large a by the 
following ansatz for an 'extended' parametrization 

r2 



and require Wi > w 2 and C 2 > . The associated profile function is 



(87) 



A E (E) 



1 + 



2C\ 
E 2 



2E'arctan ■ 



E 



Wi In 1 + 



E 2 



wi 



4C 2 ^F 




(88) 



having branch points at E = iE = iw 2 and E = iE' = iw x . In the Appendix the explicit 
expressions for the discontinuity across the cut and the weight function Pe(E) are listed. The value 
of the 'extended' profile function at E = is 



A E (0) = 1 + 2C\ 



i - c, 

w 1 



2w 



3/2 



(89) 



We determine the variational parameters C\ , to,' 



1, 2 under the constraint that the elastic threshold 



for meson-nucleon scattering is at si J u 

° thresh 



relation 



A E 



(M + m) 2 . According to Eq. (74) this requires that the 

Tfl 

= Mm H (90) 



with E = w 2 should be fulfilled during variation. 

As the imaginary part of the forward scattering amplitude below the first inelastic threshold is 
fully determined by the weight function ( see Eq. (75) ) we can now derive the total cross section in 
analytic form, in particular at threshold. After some algebra we obtain 



^totdPl = 0) = 167r 2 a 



A C\C 2 M 2 m 2 e 2m «°) 



2ttMto 



W7, 



RHw 2 ) 



(91) 
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where R(w 2 ) is the real part of the profile function at the branch point (see the Appendix). Eq. (91) 
shows that the constant total cross section at threshold is entirely due to the new term C 2 in the 
'extended' parametrization. Since we expect the overall strength C\ to be of order g 2 (see Eq. (80) ) 
the expression (91) actually is of order g 4 in the perturbative limit, in agreement with the square of the 
Born amplitude (70). This is a necessary condition for fulfilling the optical theorem and maintaining 
unitarity in our variational approach. 



5 Numerical Results and Discussion 

We first turn to the minimization of the variational functional on the pole of the propagator 

- Ml < (A 2 - 2A) M 2 + 2 (Q + V) (92) 

for the 'extended' parametrization (87). This minimization is done with respect to the 'velocity' 
parameter A and the variational parameters which enter the profile function A(E). The explicit 
expressions for the quantities and V, which are functionals of the profile function and the pseudo- 
time, can be found in (II), where also the numerical procedures are described. Again we have chosen 
M = 939 MeV and m = 140 MeV for the masses. For the present calculation the threshold constraint 
(90) was used to eliminate the parameter C 2 via Eq. (89) so that a 4-parameter minimization of Eq. 
(92) had to be performed. 

Table 2 gives the parameters of the 'extended' parametrization obtained in this way as well as some 
quantities of interest derived from them. In view of the expressions (80) and (84) for the variational 
retardation function we write the strength parameters as 

c ' = * 5& = * £ M ' < 931 
a < = *' m VfW (941 

and list the dimensionless numbers Xi,x 2 . Whereas x x stays close to its perturbative value one, the 
strength parameter x 2 for the asymptotic term in the retardation function turns out to be much larger 
even at small coupling and changes rapidly when the coupling constant a is increased. 

Comparing the numerical values of the intermediate mass Mi with the ones obtained in (II) we see 
that the 'extended' parametrization is slightly better than the 'improved' one but, of course, inferior 
to the (unconstrained) variational solution. This is more visible in Fig. 7 where 

AM? = M? - M? . . (95) 

improved 

is plotted as a function of the coupling constant a . For comparison we also have included the 
perturbative result 



M{ 



a 2 r 1 ( M 2 u 2 
M 2 - ^— / du In 1 + — ) = M 2 ( 1 - 1.0214 a ) . (96) 



pert 47r 



o 



As expected the perturbative result and the minimal value with Feynman's parametrization always lie 
above the value obtained with the 'improved' parametrization. Note that at larger coupling constants 
one has to scale these results by several orders of magnitude in order to display them in the graph. 
In contrast, the gain achieved with the 'extended' parametrization and the variational calculation is 
rather modest. As mentioned before this is due to the sensitivity of the self-energy to small values of 
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Table 2: Variational calculation for the nucleon self-energy in the Wick-Cutkosky model using the 
'extended' parametrization (87) for the retardation function. The parameters a;,, to,' , i = 1, 2 obtained 
from minimizing Eq. (92) are given as well as A for different values of the coupling constant a. They are 
constrained such that the correct elastic threshold is obtained. The dimensionless strength parameters 
Xi and x 2 are defined in Eqs. (93) and ( 94). The lower part of the table lists the mass Mi, the value 
of the profile function at E = 0, the first order residue (see (II)) and the root-mean-square radius (57) 
of the dressed particle. 





a = 0.1 


a = 0.2 


a = 0.3 


a = 0.4 


a = 0.5 


a = 0.6 


a = 0.7 


a = 0.8 


X i 


1.0093 


1.0182 


1.0296 


1.0426 


1.0600 


1.0863 


1.1266 


1.1905 




2.189 


2.564 


2.936 


3.682 


4.803 


6.363 


7.246 


4.550 


y /w 1 [MeV] 


632.7 


609.0 


584.8 


554.0 


518.0 


476.6 


437.5 


399.5 


[MeV] 


373.0 


370.0 


366.5 


362.5 


357.8 


351.8 


343.1 


324.8 


A 


0.97297 


0.94389 


0.91223 


0.87718 


0.83739 


0.79033 


0.72987 


0.62457 


M 1 [MeV] 


890.25 


839.78 


787.43 


732.97 


676.20 


616.98 


555.57 


493.45 


A(0) 


1.0150 


1.0321 


1.0518 


1.0751 


1.1037 


1.1415 


1.2002 


1.3393 


zw 


0.96088 


0.91919 


0.87429 


0.82523 


0.77042 


0.70679 


0.62687 


0.49289 


(r 2 ) 1/2 [fm] 


0.017 


0.025 


0.032 


0.038 


0.043 


0.052 


0.076 


0.135 



a so that the completely different asymptotic behaviour of the retardation function which is built into 
the 'extended' parametrization is not fully reflected in the value of M\. This is, of course, well known 
from applications of the Ritz variational principle in quantum mechanics : even very refined wave 
functions lower the ground state energy only by a small amount compared to crude ones. However, 
calculation of other observables may lead to very different results. In the present case this phenomenon 
is amplified by the need of an analytic continuation to Minkowski space in which scattering takes place. 

With the variational parameters of the 'extended' calculation fixed we can now calculate the 
imaginary part of the scattering amplitude very easily from Eq. (75). The first column of Table 3 
gives the total cross section (72) at p = obtained from the optical theorem for various coupling 
constants and Fig. 8 displays f"tot ( I P I ) f° r ^ wo different couplings below the first inelastic threshold. 
It should be noted that the threshold condition (74) associated with the first branch point E = w 2 
automatically gives the correct higher (inelastic) thresholds. However, there are additional higher 
thresholds coming from the second branch point E' = Wi in the profile function of the 'extended' 
parametrization. Since Wi is decreasing for larger coupling constants (see Table 2) these additional 
thresholds may even come to lie below the first inelastic threshold ( at |p| = 214.2 MeV ) which is 
seen, e.g., as a cusp in the total cross section for a = 0.6 near |p| = 180 MeV. It is, of course, possible 
to fix these additional thresholds at the correct physical values in a similar way as was done for the 
elastic threshold. We will not do this in the present work but concentrate on the kinematical region 
close to the elastic threshold. 

It is much more demanding to calculate the real part of the scattering amplitude at threshold in 
a reliable way. Numerical problems do not arise, of course, in the crossed amplitude which has the 
euclidean proper time integral representation (61) but in the evaluation of the oscillatory integrals 
of the direct amplitude (67). The expansion (65) is of no help since all powers of p(E) contribute 
to the real part in a given interval between thresholds. In addition, one would have to perform high 
dimensional principal value integrals numerically which is not a very promising procedure. Rather, we 
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Figure 7: Minimum of the variational functional (92) relative to the 'improved' parametrization as 
a function of the coupling constant for various parametrizations. The notation is as in Fig. 6. In 
addition, 'PTV denotes the result from first-order perturbation theory and 'ext' the one from the 
'extended' parameterization (87). 



0.5 




Figure 8: Total cross section for two coupling constants as a function of the center-of-mass momentum 
|p| . The values for a = 0.2 have been multiplied by a factor 10. The 'extended' parametrization of 
the retardation function has been employed. 



21 



have employed the subroutines DQDAWF from the IMSL Mathematical Library and D01ASF from 
the NAG Fortran Library to perform the required sine and cosine transforms (over an infinite range) 
numerically. The latter routine was used to calculate £(ir) with high accuracy whereas the former one 
served to evaluate the r-integral in (67) with less accuracy. To obtain a properly convergent integral 
the corresponding Born term first had to be subtracted, i.e. the decomposition 

-( 2m 2 -t)C(ir) _ -i , I" -(2m 2 -t)£(ir 



1+ [e-^ 2m ~ t >^ T > - 1 j (97) 

was made in the integrand. While the first term generates the Born term for the direct diagram the 
last one is now amenable to numerical integration as £(ir) goes to zero for large values of r (see Eq. 
(A. 27) ). We have checked the numerical stability of our program by treating the crossed diagram in 
the same fashion and comparing the results with the euclidean proper time representation (61). Table 
3 also contains the values for 



CT eldPl = °) = 7T~ \ Re T(s 1 t = 0)] 2 (98) 

167TS L J 



at threshold where s = (M + to) 2 and the imaginary part of the scattering amplitude vanishes. For 
comparison we also list the corresponding Born cross section (see Eq. (70) ) 



<r™ b °* B °™ (|p| = 0) = r- = 0.4226a 2 [fm 2 ]. (99) 

(M+to) 2 (1-^ 



Table 3: Total and elastic meson-nucleon cross sections at threshold in the 'extended' parametrization 
for different values of the coupling constant a. The total cross section (91) has been evaluated 
assuming the optical theorem whereas the elastic cross section has been obtained from the square of 
the scattering amplitude at threshold. For comparison the elastic Born cross section (99) is also given. 
The last column lists the 'unitarity ratio' between elastic and total cross section. 



a 


CT tot [fm 2 ] 


a el [fm 2 ] 


^Born [fm 2] 


o'elMot 


0.1 


0.0024 


0.0045 


0.0042 


1.889 


0.2 


0.0119 


0.0195 


0.0169 


1.634 


0.3 


0.0330 


0.0478 


0.0380 


1.448 


0.4 


0.0802 


0.0942 


0.0676 


1.176 


0.5 


0.183 


0.168 


0.106 


0.916 


0.6 


0.409 


0.288 


0.152 


0.704 


0.65 


0.585 


0.375 


0.179 


0.640 


0.7 


0.779 


0.495 


0.207 


0.636 


0.75 


0.929 


0.665 


0.238 


0.716 


0.8 


0.835 


0.956 


0.270 


1.145 



We observe that at larger coupling constant the Born cross section is enhanced by up to a factor 3 
due to final state interactions, vertex corrections and self-energy effects which are all (approximately) 
contained in our result. Most important are the vertex corrections because the enhancement can be 
explained, with an accuracy of better than 10 % , by replacing the coupling constant g which enters 
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Eq. (99) by the effective coupling constant </ e ff. More precisely but numerically nearly identical is 
the replacement 2g — > G 2 \(q 2 = — m 2 ) (see Eq. (54)) which leads to 

-el dPl = 0) - 4 ™) At ma , . (100 ) 



Since the momentum transfer is small compared to the nucleon mass the form factor stays very close 
to one and the main enhancement effect comes from the factor , i.e. the effective coupling constant. 
The 'unitarity ratio' 

v. = '*M = Z (101) 

^totdPl = 0) 

decreases from 1.9 to 0.6 when the coupling is varied from a = 0.1 to a = 0.7 before rising again. 
These values can be understood semi-quantitatively in the following way: we first rewrite Eq. (91) as 

a tot (|p| = 0) = na A (M + f)4 e (102) 

where we have used the definitions (93, 94) for the strength parameters. We then assume that the 
elastic cross section can be approximated by Eq. (100) which gives for the ratio 

V,«±. (1 tr )4 ^, • • S# . (103) 



y 1 ^ m) V 4M2 / 

Note that the large enhancement due to fourth power of A has cancelled. In addition, the second 
factor is l + 0{m 2 /M 2 ) , the form factor is practically one and the last factor also turns out to be very 
close to one (except at a = 0.8 where it is 1.17 ). Thus for nearly all accessible coupling constants 
one has the simple result 

U « — . (104) 

X 2 

From Table 2 we see that the dimensionless parameter x 2 grows from 2.2 at small coupling to over 
7 at a = 0.7 before declining again and that the approximation (104) accounts rather well for the 
values c e }/ctot listed in Table 3. 

The reverse procedure also works satisfactorily as can be seen in Fig. 9 : Here we have plotted 
the unitarity ratio for the 'extended' parametrization together with the results from a variational 
calculation in which x 2 has been fixed to the value x 2 = 4. (This leads to a minimal value of 
the variational functional which is nearly as good as the one from the unconstrained 'extended' 
parametrization.) Except for coupling constants close to the critical coupling we now observe equality 
of elastic and total cross section to a much better degree. It is clear that a fine tuning of the parameter 
x 2 could lead to a completely unitary result, at least at |p| = . 

However, our aim is not to unitarize the scattering amplitude but to evaluate the imaginary part of 
the amplitude as a prediction of our variational approach. To what extent unitarity is fulfilled 3 thus 
serves as a severe test of our approximation scheme, in particular near threshold. In contrast, imposing 
unitarity, e.g. by considering the Born terms as K-matrix elements [16], is an ad-hoc procedure which 
is applied on top of an approximate calculation which violates unitarity to a much larger extent. 

We also have performed a partial-wave projection of the scattering amplitude. Given the particular 
t-dependence of Eqs. (61) and (67) this can be done analytically. Fig. 10 shows the Argand diagram 
for the s-wave at a = 0.5. We observe that up to |p| ~ 60 MeV/c the s-wave amplitude remains on 
the unitarity circle before appreciable deviations occur. 



3 Of course, here we disregard the instabiltity of the Wick-Cutkosky model and assume that the theory is unitary 
below the critical coupling. 
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a 



Figure 9: Ratio of elastic to total cross section at threshold as a function of the coupling constant a. 
The full line gives the result using the 'extended' parametrization from Table 2, whereas the dotted 
line follows from a variational calculation in which the strength parameter x 2 = 4 has been kept fixed. 
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Figure 10: Argand plot for the s-wave scattering amplitude with the 'extended' parametrization of the 
retardation function at a = 0.5. The triangles denote the values for different center-of-mass momenta 
from to 140 MeV/c in intervals of 10 MeV/c. 
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6 Summary and Conclusion 



We have extended the polaron variational approach in the scalar Wick-Cutkosky model to the scatter- 
ing or absorption of an arbitrary number of mesons from the dressed nucleon. The retarded quadratic 
trial action whose parameters have been fixed variationally on the pole of the nucleon propagator was 
employed instead of the exact effective action. This constitutes the zeroth order approximation in a 
systematic expansion of the Green functions around the trial action and is similar to a calculation 
in quantum mechanics when the variational wavefunction determined from minimizing the energy is 
used to calculate other observables. 

We have seen that already this lowest order gives sensible results: although only agreement with the 
tree level calculations is assured in the perturbative limit, the variation of the parameters effectively 
sums up parts of higher diagrams up to all orders. A nice example of this is the (2 + ra)-point function 
in Section 3. This expression contains pieces which may be identified with Feynman diagrams of 
arbitrary complexity and to any order in the coupling. A look at the explicit expressions also shows 
that the zeroth order variational approximation exponentiates lowest order results in a particular way. 
This is a welcome feature since exponentiation is a frequently used recipe to extend the range of 
validity of perturbation theory. As a caveat one should add, however, that the off-shell behaviour of 
the zeroth order propagator was found to be unsatisfactory if the on-shell variational parameters are 
used to extrapolate away from the nucleon pole. No problems arise if only truncated (on-shell) Green 
functions are considered as done in this work. 

Independent of the (not very realistic) model field theory which we consider here the general 
expressions even may be used in phenomenological applications by parametrizing the functions £(r) 
or p(E) which fully determine all on-shell Green functions. Once this is done for elastic meson-nucleon 
scattering all other multi-meson processes could be predicted. 

In the present work we first have studied in some detail the vertex function for the absorption 
of a virtual meson on the dressed nucleon. Given that the trial action is quadratic in the nucleon 
trajectories it is not surprising that the corresponding form factor turned out to be gaussian. For the 
radius of the dressed particle we obtained a similar expression as in the polaron case. Since there is no 
tree-level radius the numerical results showed some differences between the various parametrizations 
which enter the trial action. 

We then concentrated on the zeroth order 'Compton' amplitude for meson-nucleon elastic scatter- 
ing which has a much richer physical content. This required an analytic continuation of the variational 
results obtained in euclidean space back into Minkowski space. We have shown that the key to a suc- 
cessful description of scattering at threshold is the proper form of the retardation function f(cr) which 
multiplies the quadratic trial action. For example, Feynman's classic parametrization used for the 
polaron, also employed by Mano [8] for the nucleon self-energy, is ruled out as it gives rise to a totally 
inappropriate analytic structure in the complex energy plane. Incidentally Mano himself writes in the 
conclusion of his work: "We also note that this method can be extended, though the accuracy of the 
result may be not very high, to another problem such as the scattering of the meson by the nucleon 
by using the best estimate of the real action." 

We already had found in previous work (I, II) that a good variational calculation of the self- 
energy requires an l/<7 2 -singularity for small a. We now find in addition that scattering near the 
elastic threshold demands a specific behaviour of the retardation function for asymptotic values of the 
proper time a. Remarkably the variational solution for the retardation function which was derived on 
the nucleon pole already contains that information (see Eq. (84) ) and has guided us to the appropriate 
form of the retardation function in both limits. We have incorporated the small- and large-u behaviour 
in an 'extended' parametrization that gives the correct elastic threshold and an imaginary part of the 
scattering amplitude which grows linearly with the center-of-mass momentum away from threshold. 
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Such a momentum-dependence leads to the expected constant total cross section at threshold. In 
addition, the 'extended' parametrization gives a lower value of the variational functional than the 
previously studied parametrizations (see Fig. 7). As our variational principle is a minimum principle, 
this minimum value is a clear measure of the quality of the corresponding ansatz. 

By means of the optical theorem we have calculated the total cross section and compared it with 
the integrated elastic cross section. The latter shows a considerable enhancement over the Born cross 
section at larger coupling. This is mainly due to vertex corrections which give rise to a larger effective 
coupling constant. At threshold we found a ratio of elastic to total cross section between 1.9 and 0.6 
depending on the coupling constant. By a simple analytic approximation we were able to show that 
this ratio is mainly determined by the strength of the asymptotic part of the retardation function 
and that a slight readjustement gives nearly unitary results. It is remarkable that threshold position, 
threshold behaviour of the total cross section and unitarity basically can fix all quantities in the 
asymptotic form (85) of the retardation function. 

Although unitarity requires strict equality of elastic and total cross section below the first inelastic 
threshold (if the instability of the model is disregarded), we still consider the numerical result satis- 
factory in several respects. First, it is a prediction of our zeroth order variational principle without 
invoking any unitarizing procedure. Second, the 'extended' parametrization is still not the optimal 
variational solution as Fig. 7 shows. Given the sensitivity of the analytic continuation procedure to 
small changes in the retardation function one may expect a further improvement of the unitarity ratio 
when more refined ansatze are used. A solution of the variational equations with the correct elastic 
threshold as constraint would be the optimal procedure if the analytic continuation into the complex 
plane could subsequently be performed. 

However, a more promising strategy is to extend the variational principle to the (2 + ra)-point 
function itself. This requires the consistent amputation of the dressed nucleon propagators and auto- 
matically leads to agreement with first order perturbation theory for small coupling constants. In a 
future publication we will show that such an extension is indeed possible at least for the simple model 
field theory which we have considered up to now. Further corrections in powers of the difference 
between the exact effective action and the trial action can then be calculated in a similar way as for 
the polaron problem [17]. The variational principle in the particle representation of field theory thus 
leads to a systematic sequence of nonperturbative approximations. 
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Appendix : Weight function for £(t) 



Here we give explicit expressions for the weight function p(E) which determines the function £(r) 



via 



/>00 

£(r) = / dE p(E) e~ Er 

JEo 



(A.l) 



for the various parametrizations of the retardation function and derive some general properties of p 
and £ in terms of the associated retardation function. 

We start with Feynman's profile function. Although it does not have the analytic structure in the 
complex Emplane which we we have assumed it is nevertheless possible to apply Eq. (64). This gives 



Pf(E) 



v - w 



4v 3 



S(E-v) 



(A.2) 



and obviously yields Eq. (50) when substituted into Eq. (A.l). 

For profile functions A(E) which have a cut running along the imaginary E-axis from iE to 
infinity we introduce 



G(E) : 
R{E) : 

Then Eq. (64) for the weight function reads 



P(E) 



1 



iplm A(iE±e) 
ReA(iE±e) . 



G(E) 



2ttE 2 R 2 (E) + G 2 (E) 



and we have to study the functions G(E) and R(E). 

For the 'improved' parametrization we write the profile function (15) as 



ME) = 1 



v 2 — w 2 1 
w W 2 



, , , iE\ . , f iE 

(w - lE In 1 + (w + iE In 1 + — 

w J \ w 



and for E > we can read off the following expressions 

v 2 — w 2 E — w 



Gi{E) 
Ri{E) 



1 + 



w E 2 
V 2 — w 2 1 

w W 2 



O(E-w) 



(w + E) In 1 + 



E 



w 



(E - w) In 



E 



w 



Note that Gi(E) > and that near threshold 



G/(E) — > 7T — (E - w) 



Ri{E) 



1 + 2- 



ln2 . 



For the 'extended' parametrization we write Eq. (88) as 



A E (E) = 1 



2d 

E 2 



iE 



(w 1 - iE) In ( 1 - — ) + ( Wl + iE) In 1 + 

w 1 J \ w 1 



iE 



+ 2 C 2V ^ ( Vw 2 + iE + Vw 2 ~ iE - 2^ 



(A.3) 
(A.4) 



(A.5) 



(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A.10) 



(A.ll) 
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and obtain 



G E {E) 



Re{E) 



2nd 

E 2 



(E - Wl ) Q(E - Wl ) 



2C- 



1 + 



2C\ 
E 2 



(w 1 + E) In ( 1 + — ) - (E - t»i) In 

Wi 



2 VE - w 2 @(E - w 2 ) 



w 1 



(A.12) 



+2C 2V ^ ( V™ 2 + £ + Vw 2 - E Q(w 2 - E) - 2^/w~ 2 



In particular, 



G E (E) 



Re(E) 



E^tw 2 



- w 2 



(A.13) 



(A.14) 



1 + 



2C\ 



{w 1 + w 2 ) In ( 1 + — j _|_ ( Wl - w 2 ) In ( 1 - — 
V w 1 J \ w 1 



w 2 \ 



-2(2 - V2)C 2V f- 



7TW 2 



(A.15) 



Note that Ge(E) is negative at least between the first and the second branch point. 

Some general properties of the functions G(E), R(E), A{E) and f(cr) are worthwhile to be noticed 
and illustrated by the particular examples we have given above. The real part R(E) is, of course, 
related to G(E) by a dispersion relation 



(A.16) 



which can be further simplified by noting that G(E') is odd and therefore R(E) is even. The profile 
function itself may be expressed by G(E) as 



2 r°° 
A(E) = 1 + - / 
vr Jeo 



dE 1 E ' „ G(E') . 

E' 2 + E 2 v ' 



At small E this becomes 



o />oo 1 o />oo 1 

1 + - / dtf ^ - i? 2 - / dE' — G(E') + 

7T 7_E 7T J_E E' a 



(A.17) 



(A.18) 



If G(E') is negative near threshold (as is the case for the 'extended' parametrization) then the coeffi- 
cient multiplying the — E 2 term may become negative as well. The profile function then will rise at 
small E from the value 

A(0)=1+- / dE-G(E). (A.19) 

7T Je E 

We can also express the inverse of the profile function by the discontinuity of l/A(iE) across the cut, 
i.e. by p(E). In particular, at E = we find 



A(0) 



1 - 4 / dEEp(E) 

E 



(A.20) 



assuming that A(E) has no zeroes in the upper half-plane. This can be checked numerically by 
comparing Eq. (A.20) with either Eq. (A.19) or with the explicit analytic expression, if available. 
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For the 'extended' parametrization we have verified the equality to be better then 1 part in 10 8 for 
all coupling constants. 

It is also easy to see that the retardation function f(cr) is the Laplace transform of E 2 G(E). 
Inversion then gives a direct relation between G(E) and the retardation function 



1 /. + 00 

G{E) = — da f(ia + e)e* E \ E>0. (A.21) 

J — oo 



The large E'-behaviour of G(E) is linked to the the small u-behaviour of f(cr) and vice versa. For 
example, the l/<7 2 -singularity of the realistic retardation functions leads to the asymptotic 1/E- 
behaviour of G(E) which is observed in Eqs. (A. 7) and (A. 12). Conversely, Eq. (A.21) may be used 
to show that a retardation function which asymptotically behaves like 

f(a) £ e- E °° (A.22) 



leads to a behaviour near threshold like 

G{E) MW) {e " jB ° r_1 @{E " Eo) ' (A ' 23) 



With Eq. (A. 5) this implies that the corresponding weight function has the threshold behaviour 



where R(E ) is the real part of the profile amplitude at the first branch point. 

Finally we give the asymptotic expansions for £(r) for r being either small or large. Assuming 
the canonical <7 _2 -behaviour of the retardation function at small a one finds in the first case 

Z(r) ^ m - (l - ^) \ - 6 r 2 In r + (9(r 2 ) (A.25) 

where 

6 = I lim i a 2 f(a). (A.26) 

Z <7->0 

The logarithmic term is due to the l/E^-fall off of p(E) for large E which does not allow a naive 
expansion of the exponential in Eq. (A.l). It produces a cut for £(r) on the negative real r-axis. The 
linear term in r has been expressed by A with the help of Eq. (A. 20). This is in agreement with the 
relation (48) between £(r) and the pseudotime and the small-r-behaviour (19) of the latter. 
For t — > oo the threshold behaviour (A. 24) of the weight function is relevant and leads to 

For purely imaginary r (which is needed in the analytic continuation) this results in a relatively slow 
and oscillating decrease at infinity. 
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